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THE  DISPERSION  OF  FIRE  SUPPRESSION  AGENTS  DISCHARGED  FROM  HIGH 
PRESSURE  VESSELS:  ESTABLISHING  INITIAIVBOUNDARY  CONDITIONS  FOR  THE 

FLOW  OUTSIDE  THE  VESSEL 


Leonard  Y.  Cooper 
Fire  Modeling  Group 
Building  and  Fire  Research  Laboratory 
National  Institute  of  Standards  and  Technology 


ABSTRACT 

This  work  reports  on  part  of  an  effort  to  study  the  dispersion  and  extinguishment  effectiveness  of 
Halon  and  Halon-altemative  fire  extinguishment  agents  discharged  from  N2-pressurized  vessels.  In 
the  systems  under  consideration,  as  the  agent  exits  fi'om  the  vessel,  thermodynamic  and  fluid-dynamic 
instabilities  lead  to  flashing  and  break-up  of  the  agent  into  a two-phase  droplet/gaseous  jet  mixture. 
This  occurs  in  a transition  region  relatively  close  to  the  vessel  exit  orifice/nozzle.  Downstream  of  this 
region  the  two-phase  agent  jet  then  mixes  with  the  ambient  air  environment  and  is  dispersed  in  the 
protected  space. 

A mathematical  model  has  been  developed  previously  to  simulate  the  time-dependent  discharge  of 
the  agent  from  the  pressure  vessel.  Using  the  output  of  this  model  and  thermodynamic  and  fluid- 
dynamic  considerations  of  the  phenomena  in  the  transition  section,  the  present  work  develops  a 
method  for  determining  a set  of  initial/boundary  conditions  at  an  initial  section  of  the  jet,  downstream 
of  the  transition  region.  These  initial/boundary  conditions  are  in  a form  that  can  be  used  to 
formulate  and  solve  the  problem  of  the  development  and  dispersal  of  the  ensuing  mixed  air/two- 
phase-agent  jet. 

Example  applications  of  the  developed  methodology  are  presented.  These  are  for  agent  discharge 
fi'om  a half-liter  cylindrical  discharge  vessel  with  a circular  discharge  nozzle/orifice  of  diameter 
0.019m.  Simulations  involve  discharge  of  the  vessel  when  it  is  half-filled  with  either  Freon  22  or 
Halon  1301  and  then  pressurized  with  N2  to  41.37xlO^Pa  (600psi). 


Key  Words:  agents,  aircraft  fire  safety,  discharge,  fire  extinguishment,  fire  safety,  halon,  halon 
alternatives 


BACKGROUND  AND  PURPOSE 


The  Building  and  Fire  Research  Laboratory  (BFRL)  of  the  National  Institute  of  Standards  and 
Technology  (NIST)  is  carrying  out  a Program  for  the  US  Air  Force  to  study  the  dispersion  and 
extinguishment  effectiveness  of  Halon  and  Halon-altemative  fire  extinguishment  agents  discharged 
from  N2-pressurized  vessels.  This  is  being  carried  out  in  support  of  an  effort  to  advance  the 
technology  of  fire  safety  in  US  Air  Force  aircraft. 

One  objective  of  the  program  is  to  measure  and/or  predict  the  dispersion  of  a discharging  agent 
throughout  a protected  space.  Time-dependent  agent  concentrations  associated  with  a particular 
agent  discharge  and  threat  scenario  would  be  used  as  a basis  for  estimating  extinguishment 
effectiveness. 

This  work  focuses  on  one  aspect  of  the  problem  of  predicting  agent  dispersion  by  means  of 
mathematical/computer  modeling. 

The  strategy  to  predict  agent  dispersion  characteristics  involves  two  basic  mathematical  model 
components:  1)  a model  to  simulate  the  time-dependent  discharge  of  the  agent  from  the  high- 
pressure  discharge  vessel;  and  2)  a model  to  simulate  the  development  and  dispersal  of  the  ensuing 
mixed  air/two-phase-agent  jet. 

A critical  element  in  exercising  any  component-2  model  is  to  establish  a set  of  initial/boundary 
conditions  for  the  jet.  This  must  be  at  a location  of  the  jet  axis  near  the  exit  section  of  the  discharge 
vessel.  The  initial/boundary  conditions  would  be  derived  from  the  predictions  of  the  component- 1 
model.  The  purpose  of  this  work  is  to  develop  a method  to  estimate  the  initial  conditions  from  the 
predictions  of  the  component-1  model  described  in  [1]. 


A MODEL  OF  THE  DISCHARGE  PROCESS 

In  response  to  the  above-mentioned  Program  objectives  a component- 1 -type  mathematical  model  of 
agent  discharge  was  developed  in  [1].  The  model  is  designed  to  simulate  agent  discharge  for  test  or 
field-deployed  system  configurations  depicted  in  Figure  1.  This  involves  a cylindrical  discharge  vessel 
with  a short,  circular,  exit  nozzle/orifice.  The  vessel  is  filled  with  pure  agent  (part  liquid  and  part 
gas).  The  agent  is  then  pressurized  with  N2  to  some  prescribed  value.  Liquid  agent  and  some 
dissolved  N2  is  in  the  lower  part  of  the  vessel.  Above  this  is  a mixture  of  gaseous  agent  and  N2.  In 
the  case  of  a field-deployed  system,  liquid-agent  discharge  is  initiated  by  actuation  of  an  explosive  cap 
sealing  the  exit  nozzle/orifice. 

A Figure-1  test  configuration,  which  would  be  used  to  simulate  field-deployed  system  discharges  under 
laboratory  conditions,  is  equipped  with  the  indicated  high-pressure  N2  holding  tank.  This  is 
connected  to  the  discharge  vessel  via  an  orifice.  This  system  simulates  the  field-deployed  system  by 
using  a rupture  diaphragm  over  the  exit  nozzle/orifice  (nominal  rupture  pressure  same  as  the  pre- 
discharge pressure  of  the  field-deployed  system)  rather  than  an  explosive  cap.  An  experimental  run 
begins  with  the  onset  of  through-orifice  N2  flow  from  the  holding  tank.  The  vessel  is  pressurized  to 
the  point  of  diaphragm  rupture  and  this  is  immediately  followed  by  vessel  discharge  to  the 
atmospheric  pressure  environment  Note  that  systems  of  interest  in  the  present  study  typically  involve 
initial  discharge  vessel  pressures,  P^y,  of  several  tens  of  atmospheres,  where,  in  the  "upside  down" 
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Figure- 1 cxjnfiguration  and  at  the  time  that  liquid  discharge  is  being  completed,  P^y  is  still  large 
enough  to  preclude  flashing  of  the  liquid  in  the  vessel. 

The  reader  is  referred  to  Reference  [1]  for  a detailed  description  of  the  model  and  for  results  of 
simulations  carried  out  to  establish  an  experimental  design  and  procedure  which  would  be  expected 
to  closely  simulate  field-deployed  system  discharges. 

The  model  simulates  the  time-dependent  discharge  of  Figure-l-type  systems  up  to  the  time  that  the 
last  of  the  liquid  agent  is  ejected  from  the  vessel.  For  a specified  set  of  geometric  system  parameters, 
agent  material,  and  initial  conditions,  key  results  of  a model  simulation  are  the  rate  of  liquid  agent 
discharged  from  the  vessel  and  the  vessel  pressure.  This  work  develops  a procedure  to  use  these 
results  along  with  other  considerations  to  estimate  initial/boundary  conditions  for  a component-2-type 
model  simulation  of  the  development  and  dispersal  of  the  ensuing  mixed  air/two-phase-agent  jet. 


MODELING  THE  EARLY  DEVELOPMENT  OF  THE  JET:  FROM  THE  EXIT  SECTION  OF  THE 
DISCHARGE  VESSEL  TO  A NEARBY  SECTION  IN  A STATE  OF  THERMODYNAMIC 
EQUILIBRIUM  AND  COMPLETED  DROPLET  FORMATION 

General  Considerations 

It  is  assumed  that  within  the  discharge  vessel  and  upstream  of  the  exit  nozzle/orifice  the  velocity  of 
the  liquid  agent  is  so  small  that  its  kinetic  energy  can  be  neglected  in  considerations  of  its 
thermodynamic  state.  Therefore,  the  thermodynamic  state  upstream  of  the  nozzle/orifice  position  can 
be  estimated  by  the  pressure  in  the  discharge  vessel,  Pj)y,  and  the  temperature  there  of  the  liquid 
agent,  Tj)y^.  Note  that  in  most  applications  it  is  expected  that  Tpy^  is  well  approximated  by 
the  ambient  temperature  outside  the  vessel.  A point  representing  the  state  of  the  liquid  agent 
in  the  vessel  is  indicated  in  the  sketched  pressure-enthalpy  (P,  h)  diagram  of  Figure  2. 

Because  of  the  relatively-short  nozzle\orifice  design  under  consideration,  the  time  interval  for  the 
liquid  to  approach  it,  pass  through  it,  and  enter  the  outside  environment  is  very  small.  For  conditions 
related  to  the  NIST  Program  this  is  of  the  order  of  lO’^s.  During  this  time  interval  the  pressure  of 
the  traversing  liquid  is  reduced  from  the  high  P^y  values  inside  the  vessel,  of  the  order  of  several 
tens  of  atmospheres,  to  pressures  of  the  order  of  P^mb’  atmosphere. 

As  discharging  liquid  enters  and  traverses  the  region  of  the  nozzle/orifice,  it  will  enter  and  move 
along  super-heated  metastable  thermodynamic  states.  While  dependent  on  its  changing  thermody- 
namic state  in  the  vessel,  (P^v*  "^dval)’  particular  thermodynamic  properties,  as  the  agent 

material  penetrates  the  ambient  environment  it  approaches,  and  likely  achieves  the  P^mb  pressure 
while  still  in  its  metastable  liquid  state.  Thus,  downstream  of  a vena  contracta  the  liquid  agent 
develops  into  a near-uniform-radius  liquid  jet,  r^j,  which  can  be  described  following  traditional 
incompressible  fluid-dynamic  considerations.  See  Figure  3.  For  the  relatively  large-vapor-pressure 
agent  materials  of  interest  here,  a combination  of  fluid-dynamic  and  thermodynamic  instabilities  will 
then  lead  to  breakup  into  small  droplets  and  flashing  (i.e.,  rapid  evaporation  to  a two-phase  equilibri- 
um thermodynamic  state)  of  the  metastable  liquid  jet  [2]. 

The  metastable  liquid  is  expected  to  move  through  thermodynamic  states  on  near-isentropic  paths. 
However,  as  discussed  in  [1],  for  a given  agent  and  initial  state  conditions,  (Pdv>  "^dv,al)’  ^ 
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possible  that  the  pressure  along  such  a path  will  not  reach  prior  to  an  intersection  with  the 

agent’s  spinoidal  curve.  Moreover,  if  such  an  intersection  does  occur,  then  at  that  instant 
spontaneous  nucleation  in  the  liquid  is  to  be  expected  [3,  4].  It  is  conjectured  that  this  would  lead 
to  near-explosive  breakup  and  flashing  of  the  jet,  all  of  this  being  initiated  within  a jet  penetration 
depth  (into  the  outside  environment)  of  the  order  of  a single  nozzle/orifice  diameter. 

Whatever  the  thermodynamic-state  path  of  the  discharging  metastable  liquid,  it  is  reasonable  to 
assume  that  a relatively  short  distance  downstream  of  the  nozzle/orifice  a thermodynamic  equilibrium 
state  and  an  end  to  the  droplet-formation  phenomenon  will  be  achieved.  It  is  also  reasonable  to 
assume  that,  as  the  agent  flows  further  downstream,  developing  as  a mixed,  two-phase,  agent/air  jet, 
thermodynamic  equilibrium  will  be  maintained,  and  droplet  collision  and  agglomeration  will  not  play 
an  important  role  in  jet  dynamics. 


The  Initial  Section  of  the  Jet 

The  position  along  the  jet  axis  where  both  approximate  thermodynamic  equilibrium  and  completion 
of  droplet  formation  is  first  achieved  will  be  denoted  here  as  the  initial  position  of  the  jet.  Associated 
with  this  position  is  the  plane  section  normal  to  the  jet  axis,  called  the  initial  section  of  the  jet. 
Similarly,  the  thermodynamic  state  and  properties  of  the  jet  at  the  initial  position  will  be  denoted  as 
its  initial  state  and  initial  properties.  The  region  of  the  flow  between  the  nozzle/orifice  discharge 
section  and  the  initial  section  will  be  referred  to  as  the  transition  region  of  the  jet. 

A detailed  description  of  the  initial  state  and  properties  of  the  jet  would  be  extremely  complicated. 
This  would  include  the  variation  across  the  initial  section  of  air/gaseous-agent  concentrations, 
temperature,  and  velocity,  and  of  the  size,  velocity,  and  temperature  distributions  of  the  liquid  agent 
droplets.  Note  also  that  in  the  present  application  all  of  these  variables  of  the  initial  jet  flow  field 
would  be  time-dependent.  In  general,  an  accurate  detailed  description  of  the  initial  state  of  the  jet, 
whether  determined  by  theoretical  or  experimental  means,  is  either  impractical  or  beyond  the  current 
state  of  technology. 

Although  it  is  not  possible  to  provide  a detailed  description  of  the  initial  state  of  the  jet,  it  is 
reasonable  to  expect  that  an  achievable  approximate  description,  suitable  for  use  in  a component-2- 
type  model,  would  yield  good  approximations  of  the  dispersing-agent-flow-field  problem.  It  is  the 
purpose  of  this  work  to  provide  such  an  approximate  description. 


Assumptions  of  the  Phenomena  in  the  IVansition  Region:  From  the  Nozzle/Orifice  to  the  Initial 
Section 

The  dynamic  processes  that  initiate  breakup  and  flashing  of  the  metastable  liquid  jet  and  bring  it  to 
its  initial  state  are  not  completely  understood.  At  one  extreme,  it  is  possible  that  the  above- 
mentioned  instabilities  lead  to  violent  breakup  and  flashing  of  the  liquid  jet  immediately  upon  leaving 
the  exit  nozzle/orifice.  Such  behavior  was  observed,  for  example,  in  some  experiments  of  [5].  Also, 
toward  the  completion  of  the  liquid  discharge,  this  was  also  observed  in  recent  experiments  at  NIST 
[6]. 


4 


The  metastable  state  can  also  persist  even  as  the  liquid  jet  penetrates  relatively  deeply  into  the 
ambient  environment.  For  example,  Reference  [5]  describes  a liquid  water  jet  of  diameter  D = 
7.6(10'^)m,  initially  at  9.03xl0^Pa  (131psi)  and  415K,  which  penetrated  up  to  0.0254m,  corresponding 
to  L/D  ~ 32,  before  initiation  of  violent  breakup/flashing.  Observed  in  the  recent  NIST  experiments 
[6]  was  a liquid  Freon  22  jet,  D = 0.02m,  initially  at  41.37xlO^Pa  (600psi)  and  294K,  which  penetrated 
a distance  of  approximately  0.3m,  corresponding  to  L/D  ~ 15,  before  initiation  of  violent  break- 
up/flashing. In  this  regard  it  is  of  interest  to  note  that  in  the  literature  of  the  discharge  of  flashing 
liquids  from  high  pressure  vessels  there  is  a degree  of  conventional  wisdom  that  analyses  of  the  two- 
phase  flow  phenomena,  based  on  thermodynamic  equilibrium  states,  tend  to  yield  good  results  when 
the  nozzle/piping  from  the  inside  to  the  outside  of  a high  pressure  discharge  vessel  is  at  least  0.1m 

m- 

The  initial  section  will  be  some  incremental  distance  downstream  of  the  penetration  distance  of  the 
metastable  liquid  jet.  The  studies  of  [5]  and  the  recent  experiments  at  NIST  [6]  suggest  that  this 
incremental  distance  will  typically  not  exceed  the  order  of  the  penetration  distance  itself.  Whatever 
the  distance,  it  is  assumed  here  that  the  total  length  of  the  transition  region  is  small  enough, 
compared  to  the  characteristic  length  of  the  overall  agent  dispersal  problem  of  interest,  that  the 
processes  within  the  transition  region  can  be  treated  as  quasi-steady. 

In  the  transition  region  it  is  assumed  that  momentum  transfer  and  any  net  work  or  heat  transfer 
interactions  between  the  penetrating  agent  jet  and  the  ambient  air  environment  are  negligible.  Thus, 
consistent  with  the  analysis  of  [8]  for  steady,  two-phase,  air/evaporating-liquid  jets,  it  is  assumed  that 
the  jet  at  the  initial  section  can  be  approximated  as  consisting  only  of  a two-phase  mixture  of  agent, 
with  no  entrained  air.  Finally,  in  the  transition  section  it  is  assumed  that  increases  in  the  energy  of 
the  liquid  due  to  mechanical  work  by  surface  tension  forces  to  form  the  initial-section  droplets  are 
negligible  compared  to  the  energy  transfers  involved  in  phase-change  processes. 


The  Initial  State  - A Summary  of  Assumptions  and  Approximations  and  Their  Implications 

Regarding  an  estimate  of  the  initial  state,  the  assumptions  of  the  previous  paragraph,  some  additional 
approximations  and  assumptions,  and  their  implications  are  now  summarized. 


At  the  initial  section.  The  jet  is  defined  by  a characteristic  initial  radius  rj. 
consists  of  a pure  mixture  of  saturated-agent  vapor  and  liquid  (droplets) 
equilibrium  at  temperature,  Tj,  and  pressure,  Pj,  where 


The  region  r < rj 
in  thermodynamic 


~ ^AMB’  ^ n (1) 

~ ^SAT  “ ^SAtO* AMb)^  ^ - •'I  (2) 

and  where 
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'^SAT  ~ ^SAtC^)’  ^ - ^CR 


(3) 


is  the  saturation  temperature  of  the  agent,  a specified  function  of  pressure  P,  where  P can  not  exceed 
the  critical  pressure,  P^j^.  The  liquid-agent  component  of  the  jet  consists  of  droplets  of  uniform- 
radius  rj  j)  and  these  are  uniformly  dispersed  over  the  region  r < rj.  The  agent  droplets  and  gas  have 
the  same  uniform  axial  jet  velocity,  Vj^. 


Vi  = V,^;  r < r,  (4) 

where  Vj  = Vj(r)  is  the  velocity  distribution  of  the  jet  at  the  initial  section.  There  is  no  agent  in  r 
> rj. 

The  region  r > rj  consists  of  a uniform  ambient-temperature/pressure  environment  (P^mb*  ^amb) 
of  quiescent  air. 


“ ^AMB’ 


T,  = T 


AMB’ 


= 0;  r > r. 


(5) 


There  is  no  air  in  r < rj. 


Conservation  of  mass; 


“ dM^vAiL/dt 

(6) 

Pi  A ~ + (1  ■ 

(7) 

where  /3|  Pi^  and  Pij^Q  are  the  average  density  of  the  agent,  the  density  of  the  liquid  agent 
(droplets),  and  the  density  of  the  agent  gas  at  the  initial  section,  respectively;  Xj  is  the  quality  of  the 
two-phase  flow  at  the  initial  section,  i.e.,  the  fraction  of  a mass  of  agent  which  is  saturated  gas; 
Mdv^  is  the  instantaneous  mass  of  liquid  agent  in  the  discharge  vessel;  and  the  right  side  of  Eq. 
(6),  determined  previously  from,  e.g.,  the  discharge  model  of  [1],  is  assumed  to  be  specified. 


Conservation  of  ener2v  [and  using  Eg.  (6)1: 


^DV,AL  “ ^I,AL  *I^ALg(TsAt)  + (8) 

^ALgCTsAt)  = ■ *'l,AL  (9) 
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where  hpy^  and  hj^  are  the  specific  enthalpy  of  the  agent  liquid  in  the  discharge  vessel  and  at 
the  initial  section,  respectively;  hj^Q  is  the  specific  enthalpy  of  the  agent  gas  at  the  initial  section; 
and  h^Q  is  the  heat  of  vaporization  of  the  agent  at  the  temperature  In  Eq.  (9), 


^ALG  “ ^*alg(T)»  T ^ T^r 


(10) 


is  a specified  function  of  temperature,  T,  where  T can  not  exceed  the  critical  temperature,  T^r. 

A point  representing  the  state  of  the  agent  at  the  initial  section  is  presented  in  the  (P,  h)  diagram 
of  Figure  2.  Consistant  with  Eq.  (8),  this  reflects  the  fact  that  some  of  the  original  enthalpy  of  the 
material  in  its  high  pressure  "rest"  state  was  exchanged  for  its  now-significant  kinetic  energy. 

The  sketch  of  the  transition  region  and  of  the  initial  properties  as  presented  in  Figure  3 is  consistent 
with  the  all  the  above  assumptions/approximations. 


Droplet  size.  It  is  assumed  that  fj  ^ is  specified.  Based  on  [5]  it  is  expected  that  this  will  be  in  the 
range  10pm  - 100pm.  Photograpliic  data  during  future  NIST  experiments  should  provide  sharper 
estimates  of  this  parameter.  Also,  the  sensitivity  of  the  downstream  agent  concentrations  to  variations 
of  this  parameter  should  be  evaluated  with  component-2  model  simulations. 


Initial  jet  radius.  It  is  assumed  that  Tj  is  specified.  Based  on  photographic  data  in  [5]  and  on 
photographic  data  already  acquired  during  NIST  experiments,  rj/r^j  seems  to  be  of  the  order  of  10  - 
20.  Thus,  in  (what  is  reasonable  to  construe  visually  in  high  speed  photographs  as)  a region  of 
breakup  and  flashing  of  the  liquid  jet  to  an  initial  equilibrium  state,  the  radius  of  the  jet  increases  by 
a factor  of  10  - 20.  Again,  the  sensitivity  of  the  downstream  agent  concentrations  to  variations  of  this 
ratio  should  be  evaluated  with  component-2  model  simulations.  One  would  hope  that  such  sensitivity 
is  not  great 


Parameters  definin£  the  initial  state.  Consistent  with  the  above,  the  specified  parameters  of  the 
initial  state  are  Pj,  according  to  Eq.  (1),  rjjj,  and  Also,  for  a particular  agent  of  interest,  Tj 

is  sp)ecified  by  (2)  and  (3).  Requir^  to  complete  the  description  of  the  initial  state  are  Vj^, 
and  Xj. 

Completing  the  Description  of  the  Initial  State 

Modeling  the  liquid  agent  as  incompressible.  For  problems  of  interest  here,  agent  liquid 
temperatures  are  expected  to  be  small  enough  compared  to  Tcr  and  vary  in  a small  enough  range 
to  permit  accurate  modeling  of  the  liquid  as  being  incompressible,  with  density  p ^ and  as  having 
a constant  specific  heat,  ^AL-  Thus, 
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/^ial  “ Pal 


(11) 


The  values  of  and  will  be  taken  as 

Pal  ~ PalO^m)''  ^al  “ ^al("^al)’  = (^sat  + ^dyal)^  (^2) 

where  is  taken  to  be  the  average  temperature  of  the  liquid  in  the  discharge  vessel  and  at  the 
initial  state,  and  where  />al(T)  ^nd  C^(T)  are  specified  functions  of  T 

From  the  assumption  of  incompressibility  it  follows  that 

^DYAL  ■ " ^Al(^DY,AL  * "^SAt)  + (^DV  " PaMb)/PaL  (13) 

Modeling  the  agent  gas  as  a perfect  gas.  The  agent  gas  is  modeled  as  a perfect  gas. 

^I,AG  = I’i/C^AgTi)  = 1*AMb/(1^Ag"1^SAt)  (l*^) 

where  R^q  is  the  gas  constant 

The  solution  for  and  Xj.  As  discussed  above,  it  is  expected  that  in  most  instances  the  liquid  jet 
will  achieve  P^mb  prior  to  its  breakup/flashing.  Then,  Bernoulli’s  equation  leads  to 

~ '^u'PaaP'^dw  ’ ^ambV/^alI^^  (1^) 

where  r^j  is  determined  from  the  radius  of  the  exit  nozzle/orifice,  r^^,  and  its  orifice  discharge 
coefficient,  Cp. 

(*^lj/*’n)^  ~ (1^) 

Substituting  Eqs.  (9)  and  (13)  into  Eq.  (8)  lead  to 

~ ^Al(^DV,AL  ■ ^SAt)  (1*DV  ■ ^ambV/^al  ■ *i^alg("1^sat)  (1^) 
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Substituting  Eqs.  (7)  and  (15)  into  Eq.  (6)  leads  to 


- (rLj/ri)'*[(PDv  ■ PambV/^aiJ[(V^)(^  ' 

^ Piag^^al 

Comparing  Eqs.  (18)  and  (17)  leads  to 

+ A2X  - (Ai  + 1)  = 0 (20) 

where 

X = (xj/e)(l  - e)(rLj/ri)^;  Xj  = X£/[(l  - e)(rLj/ri)^]  (21) 

~ Pal^alC^dvm^  ■ "^sat)/(Pdv  ■ ^amb)  (22) 

A2  = [e/(l  - e)lPALhALG(TsAT)/[(PDV  ’ ^AMB)(^dhf]  (23) 

The  solution  to  Eq.  (20)  is 

X = - (A2/2)  + [(A2/2)2  + Ai  + 1]^^  (24) 


^^^lth  Eq.  (21)  and  other  previous  definitions,  Eq.  (24)  provides  the  solution  for  Xj.  This  would  be 
used,  in  turn,  to  obtain  Vj^  from  Eqs.  (17)  or  (18),  thereby  completing  the  description  of  the  initial 
state. 


Via  = 2X[(Pdv  - PambVPaJ''" 


(IT) 


A sketch  of  the  velocity  distribution  and  the  mass  flux  distributions  of  gaseous  and  liquid  agent  at  the 
initial  section  are  presented  in  Figure  4. 


ALGORITHM  FOR  DETERMINING  THE  INITIAL  STATE 
The  above  solution  can  be  implemented  with  the  following  algorithm: 
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1.  Specify  agent  and  its  thermodynamic  properties,  including  ToatCP)  and  hAi  n(T)»  according 
to  Eqs.  (3)  and  (10),  respectively,  and  Pal(T),  Cyij^(T). 

2.  Specify  rj/r^j  (expected  range  10  - 20). 

3.  Specify  P^mb’  Tamb’  and  (In  most  cases  expect  to  specify  = T^mbO 

4.  Specify  ~ ^AMB  for  r > rj.  Compute  TsatC^amb)  (^)  and  then  specify  T|  = 

Tsat  for  0 < r < r,. 

5.  Compute  from  Eq.  (10). 

6.  Compute  Tal  according  to  Eq.  (12)  and  then^AL  and  Cal  according  to  Eq.  (12)  andyOAL(T) 
and  Cal(T)»  respectively. 

7.  Compute  Pi^Q  according  to  Eq.  (14)  and  then  e according  to  Eq.  (19). 

8.  Use  a component-1  model,  e.g.,  that  of  [1],  to  find  the  time-dependent  values  of  P^y  for  a 
particular  Figure- 1 -configuration  of  interest. 

9.  For  a particular  value  of  P^jy  compute  and  A2  from  Eqs.  (22)  and  (23),  and  then  X from 
Eq.  (24). 

10.  Find  Xj  from  Eq.  (21)  and  Vj^  from  Eq.  (17’). 

EXAMPLE  APPLICATION  OF  THE  ALGORITHM 

Conditions  for  Example  Discharges  and  Temperature-Dependent  Agent  Property  Functions 

The  algorithm  was  applied  in  four  example  discharge  simulations.  These  involve  the  agents  Freon 
22  (CHCIF2)  and  Halon  1301  (CBrFj)  and  two  different  rj  values  corresponding  to  rj  = lOrjj  and 
20rLj.  The  examples  consider  time-dependent  agent  discharges  fi-om  a Figure- 1 configuration  with 
no  holding  tank.  The  geometric  parameters,  orifice  coefficient,  and  initial  conditions  assumed  in  the 
simulated  example  discharge  calculations  are  presented  in  Table  1.  These  are  the  same  parameters 
used  in  the  example  calculations  of  [1].  During  the  discharge,  the  time-dependent  value  of  P^y  was 
computed  using  the  component-l-type  model  of  [1]. 

The  material  property  functions  for  Freon  22  (CHCIF2)  and  Halon  1301  (CBrF3)  which  were  used 
in  the  example  calculations  are  from  [9].  These  are  of  the  form 


Molecular  weight,  M:  M kg  of  agent  = 1 kg-mole  of  agent  (25) 

~ exp[A  + B/(T5Ar/R)  + Cln(TsAr/R)  + "^1  - ^sat  - ^2 

(26) 
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hALG(TsAT)/[J/(kg-mole.K)]  = Mh^G(TsAT)/[J/(kg*K)] 

= A(1  - Ts^-j-zTcr)'®;  Tj  < < T2 


(27) 


/^AL(T)/(kg-mole/m^)  = (l/M)^^(T)/(kg/m^) 

= A/B^;  F = {1  + [1  - (T/R)/C]^};  < T < T2 


(28) 


C^(T)/[J/(kg-mole*K)] 


= MC^(T)/[J/(kg.K)] 

= A + B(T/R)  + C(T/R)2  + 


(29) 

D(T/R)^  Ti  < T < T2 


Complete  descriptions  of  these  functions  for  the  two  agents  are  presented  in  the  APPENDIX. 

was  obtained  as  the  root  of  Eq.  (26)  where,  according  to  Eq.  (2)  or  item  4 of  the  algorithm,  the 
left  hand  side  was  specified  as  Py^j^.  Then  was  calculated  using  Eq.  (12).  Finally,  h^Q(Ts^j), 
/?al(Tal).  and  C^(T^  were  calculated  from  Eqs.  (27),  (28),  and  (29),  respectively.  The  results 
are  presented  in  Table  2. 


Results  of  the  Example  Calculations 

The  results  of  the  example  calculations  are  tabulated  in  Table  3 (Freon  22,  rj  = lOr^j  = 

= 0.07378m),  Table  4 (Halon  1301,  rj  = lOr^j  = lOCp^^r^^  = 0.07378m),  Table  5 (Freon  22,  rj  = 
20rLj  = 20Cd^^n  = 0.1476m),  and  Thble  6 (Halon  1301,  rj  = 20ru  = 20CD^^rN  = 0.1476m. 
Included  in  the  tables  are  computed  values  of  the  variable  which  is  the  dimensionless  mass  of  liquid 
agent,  or  dimensionless  elevation  of  the  liquid  agent  in  the  discharge  vessel,  Z. 


C “ ^dv,al/^dv,al(^  — 0)  — Z/Z(t  — 0) 


(30) 


The  tabulated  results  of  Table  3 for  the  discharge  of  Freon  22  (rj  = lOr^j  = 0.07378m)  are  plotted 
in  Figures  5-8.  These  include  plots  of  the  time-dependent  values  of  Pj^y  (Figure  5),  Vj^  (Figure 
6),  ^ and  Xj  (Figure  7),  and  dMj^/dt  and  dM|^Q/dt  (Figure  8). 

A study  of  the  results  leads  to  the  following  observations: 

1.  For  both  Freon  22  and  Halon  1301,  the  times  to  complete  the  discharge  of  liquid  agent  from 
the  pressure  vessel  are  almost  identical  at  0.024  ± 0.001s  (Tables  3 - 6). 

Z At  the  initial  section,  the  quality  of  the  two-phase  agent  flow,  X|,  is  significantly  different  for 
the  two  agents  considered,  but  for  a given  agent  Xj  is  relatively  independent  of  t and  rj,  for 
the  two  initial-section  radii  values  considered  (rj  = 0.07378m  and  0.1476m).  Thus,  for  the 
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cases  considered,  Xj  = 0.31  ± 0.01  for  Freon  22  (Tables  3 and  5)  and  Xj  = 0.54  ± 0.02  for 
Halon  1301  (Tables  4 and  6). 

3.  At  the  initial  section,  for  a given  agent,  and  for  the  two  initial-section  radii  values  considered, 
the  mass  flow  rates  of  liquid  or  gaseous  agent  are  almost  identical  at  any  particular  time 
during  the  discharge.  [For  example:  for  Freon  22,  approximately  half-way  through  the 
discharge  at  t = 0.01252s,  values  of  dM|^/dt  for  rj  = 0.07378m  and  rj  = 0.1476m  are 
8.59kg/s  (Table  3)  and  8.49kg/s  (Table  5),  respectively;  and  for  Halon  1301,  approximately 
half-way  through  the  discharge  at  t = 0.01222s,  values  for  dMj^/dt  for  rj  = 0.07378m  and 
rj  = 0.1476m  are  6.95kg/s  (Table  4)  and  8.59kg/s  (Table  6),  respectively.]  However,  over  the 
course  of  the  discharge,  these  mass  flow  rates  are  reduced  to  approximately  two  thirds  of  their 
initial  value.  [For  example:  for  Freon  22  and  rj  = 0.07378,  the  values  of  dMj^/dt  at  t = 0 
and  t = 0.232fe  (near  the  end  of  the  liquid  discharge)  are  11.78kg/s  and  7.34k^s,  respectively 
(Table  3 and  Figure  5);  and  for  Halon  1301,  and  rj  = 0.07378,  the  values  of  dMj^/dt  at  t 
= 0 and  t = 0.2443s  (near  the  end  of  the  liquid  discharge)  are  9.26kg/s  and  5.88kg/s, 
respectively  (Table  4).] 
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NOMENCLATURE 


A,B,C,D 

constants,  Eqs.  (A2)  - (A5) 

^AL 

constant  specific  heat  of  liquid  agent 

orifice  coefficient 

D 

diameter  of  a liquid  jet 

F 

Eqs.  (28)  and  (A4) 

h 

specific  enthalpy 

^ALG 

specific  heat  of  vaporization  of  the  agent 

^dval 

h of  agent  liquid  in  the  discharge  vessel 

^iag 

h of  agent  gas  (liquid)  at  the  initial  section 

L 

length  of  liquid  jet 

M 

molecular  weight 

dMi^o/dt  (dMj^/dt)  mass  flow  rate  of  agent  gas  (liquid)  at  the  initial  section 
^DVAl  liquid  agent  in  the  discharge  vessel 


pressure 

^AMB 

p of  the  ambient  environment 

^CR 

critical  pressure 

^DV 

p in  the  discharge  vessel 

Pi 

p at  the  initial  section 

^AG 

gas  constant  for  the  agent  gas 

r 

radius 

r of  two-phase  agent  jet  at  the  initial  section 

•■iJD 

r of  agent  droplets  at  initial  section 

ru 

r of  the  liquid  jet  outside  the  discharge  vessel 
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^■n 

r of  the  discharge  vessel’s  exit  nozzle/orifice 

T 

absolute  temperature 

Tal 

T of  liquid  agent 

^AMB 

T of  ambient  environment 

^DVAL 

T of  liquid  agent  in  the  discharge  vessel 

Ti 

T at  the  initial  section 

"^SAT 

saturation  temperature  of  the  agent 

Ti,T2 

constant  temperatures,  Eqs.  (A2)  - (A5) 

t 

time 

^BURST 

t when  the  diaphragm  bursts 

Vi 

velocity  at  the  initial  section 

Via 

Vj  of  agent 

Vi,aIR 

Vj  of  air 

X 

Eq.  (21) 

quality  of  the  two-phase  agent  flow  at  the  initial  section 

Z 

distance  of  the  liquid/gas  agent  interface  above  the  exit  nozzle/orifice 

z 

distance  along  the  jet  axis  from  the  exit  orifice/nozzle 

Zl 

z at  the  initial  section 

e 

Eq.  (19) 

A^,  A2 

Eqs.  (21)  and  (22) 

PlA 

average  density  of  the  agent  jet  at  the  inital  section 

Piag  (Pial) 

density  of  the  agent  liquid  (gas)  at  the  initial  section 

Pal 

density  of  liquid  agent 

1 

dimensionless  elevation  of  the  liquid/gas  interface,  Eq.  (30) 
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APPENDIX:  MATERIAL  PROPERTY  FUNCTIONS  FOR  FREON  22  (CHCIF2)  AND  HALON 

1301  (CBrF3) 


The  following  material  property  functions  for  Freon  22  (CHCIF2)  and  Halon  1301  (CBrF3)  were  used 
in  the  example  calculations  [9]: 


Molecular  weight  = M kg  of  agent  = kg-mole  of  agent 


(Al) 


M 


Freon  22  fCHClFO  Halon  1301  fCBrF^^ 

86.468  148.910 


~ exp[A  + B/(Ts^,^p/R)  + Cln(T5^j/R)  + D(Ts^/R)^];  Tj  < < T2 

(A2) 


A 

B 

C 

D 

Ti 

T2 


Freon  22  (CHCIF2I 

9.7882(10) 

- 4.7735(10-3) 

- 1.2317(10) 
2.4173(10-^) 
115.73K 
369.30K 


HalgnJiOLXCBrFj} 

6.2780(10) 

- 3.3770(103) 

- 6.7461 
1.4407(10-^) 

105.15K 
340. 15K 
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- Mh^Q(T5^j)/(J/kg) 

= A(1  - Ts^j/T^r)®;  Tj  < Tg^j  < T2 


(A3) 


A 

B 

^CR  “ 
Ti  = 
T2  = 


Freon  22  rCHClF.) 

2.9617(10*^) 

3.8600(10-1) 

369.30K 

115.73K 

369.30K 


Halon  1301  rCBrF.l 

2.5070(10^) 

3.5700(10-1) 

340.15K 
105.15K 
340. 15K 


PAL(T)/(kg-mole/m^)  = p^(T)/(Mkg/m^) 

= A/B^;  F = {1  + [1  - (T/R)/C]l^};  Tj  < T < T2 


(A4) 


A 

B 

C 

D 

Ti 

T2 


R:eon22XCHClF2} 

1.5956 

2.6605(10-1) 

3.6930(102) 

2.8123(10-1) 

115.73K 

369.30K 


Halon_1301_CCBrF3} 

1.3435 

2.6995(10-1) 

3.4015(102) 

2.8006(10-1) 

105.15K 

340.15K 
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= MC^(T)/[J/(kg-K)]  (A5) 

= A + B(T/R)  + C(T/R)2  + D(T/R)^  < T < T2 


CAL(T)/[J/(kg-mole-K)] 


A 

B 

C 

D 


Freon  22  fCHClF.') 

9.4660(10^) 

3.9300(10^) 

- 8.1740(10-1) 
2.7910(10*^) 
115.73K 
369.30K 


Halon  1301  ICBrF.l 

1.1421(10^) 

- 1.3150(102) 
6.1200(10-1) 

0 

193.15K 
298. 15K 
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Thble  1. 


Geometric  parameters,  orifice  coefficient,  and  initial  conditions  for  example 
discharges  from  a Figure-l-type  configuration  with  no  holding  tank. 


Length  of  discharge  vessel  = Z^y  = 

Cross-section  area  of  discharge  vessel  = A^y  = 

Volume  of  discharge  vessel  = V^y  = 

Area  of  exit  nozzle/orifice  = Aq  = 

Radius  of  exit  nozzle/orifice  = Tq  = 

Discharge  coefficient  of  exit  nozzle/orifice  = Cj^  = 

Initial  volume  of  agent  liquid  in  discharge  vessel  = Mj^y^  j = 
Initial  temperature  of  gas  in  the  discharge  vessel  = T^yj  = 
Initial  pressure  in  discharge  vessel  = P^yi  = 

Temperature  of  agent  liquid  in  discharge  vessel  = ^^DyAL  “ 


0.2546m 

0.1963xl0-“m2 

0.500xl0-^m^ 

0.2850x1 0'^m^ 

0.9525xl0-2m 

0.6 

0.250xl0-^m^ 

294K 

41.36xlO^Pa  (600psi) 
294K 


Table  2.  Constant  properties  of  the  initial  state  for  the  example  calculations. 


Freon  22  rCHClFol 

Halon  1301  rCBrF,l 

"^SAT 

232K 

215K 

Tal 

263K 

254.5K 

^ALg(Tsat) 

0.234xl0^J/kg 

0.118xl0^J/kg 

/^al("^al) 

0.132xl0^kg/m^ 

0.180xl0^kg/m^ 

^alC^al) 

0.115xl0^J/(kg*K) 

0.809x1 0^J/(kg-K) 
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Table  3 


Initial  state  of  a Freon  22  (CHCIF2)  jet  from  Figure-l-type  vessel  with  characteris- 
tics of  Table  1,  constant  initial  state  properties  of  Table  2,  and  with  rj  = lOr^j  = 
lOCo^^rN  = 0.07378m. 


/^AL  = 1320kg/m^ 


PW  = 3.58kg/m^ 
V = 232K 


t 

[S] 

^*DV 

[Pa] 

V 

[m/s] 

*1 

[kg/s] 

[kg/s] 

0 

0.4135E+07 

1.0000 

122.73 

03022 

5.10 

11.78 

0.00179 

03663E+07 

0.9031 

115.46 

03025 

4.80 

11.07 

0.00358 

0J299E+07 

0.8117 

109.49 

03028 

4.55 

10.48 

0.00537 

0-3008E+07 

0.7248 

104.46 

03030 

434 

9.99 

0.00716 

0.2769E+07 

0.6417 

100.14 

03032 

4.16 

9.57 

0.00895 

0.2570E+07 

0.5620 

9637 

03033 

4.01 

9.20 

0.01073 

0.2401E+07 

0.4851 

93.05 

03034 

3.87 

8.88 

0.01252 

0.2255E+07 

0.4108 

90.08 

03035 

3.75 

8.59 

0.01431 

0.2128E+07 

0.3389 

87.41 

03036 

3.63 

833 

0.01610 

0.2016E+07 

0.2690 

84.99 

03037 

3.53 

8.10 

0.01789 

0.1917E+07 

0.2010 

82.78 

03038 

3.44 

7.89 

0.01968 

0.1829E+07 

0.1347 

80.76 

03039 

336 

7.69 

0.02147 

0.1749E+07 

0.0701 

78.88 

03039 

3.28 

7.51 

0.02326 

0.1677E+07 

0.0069 

77.15 

03040 

3.21 

734 
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Table  4. 


Initial  state  of  a Halon  1301  (CBrFj)  jet  from  Figure-l-type  vessel  with  characteris- 
tics of  Table  1,  constant  initial  state  properties  of  Table  2,  and  with  rj  = lOr^j  = 
lOCD^^r^  = 0.07378in. 


Pal  ~ 1800kg/m^ 


Pw  = 6.17kg/in3 
Tu  = 215K 


t 

[S] 

^DV 

[Pa] 

c 

V 

[m/s] 

[kg/s] 

dM,,^/dt 

[kg/s] 

0 

0.4135E+07 

1.0000 

142.83 

0.5181 

9.95 

9.26 

0.00204 

03687E+07 

0.9030 

135.28 

0.5205 

9.43 

8.69 

0.00407 

0.3338E+07 

0.8111 

128.99 

0.5224 

8.99 

8.22 

0.00611 

03057E+07 

0.7236 

123.64 

0.5239 

8.62 

7.83 

0.00814 

0.2826E+07 

0.6397 

118.99 

0.5252 

8.29 

7.50 

0.01018 

0.263 lE+07 

0.5590 

114.91 

0.5263 

8.01 

7.21 

0.01222 

0.2465E+07 

0.4812 

111.28 

0.5273 

7.76 

6.95 

0.01425 

0.2322E+07 

0.4058 

108.03 

0.5281 

7.53 

6.73 

0.01629 

0.2197E+07 

0-3327 

105.08 

0.5289 

732 

6.52 

0.01832 

0.2086E+07 

0.2616 

102.40 

0.5295 

7.14 

634 

0.02036 

0.1988E+07 

0.1924 

99.94 

0.5301 

7.97 

6.17 

0.02240 

0.1899E+07 

0.1248 

97.67 

0.5307 

6.81 

6.02 

0.02443 

0.1820E+07 

0.0588 

95.58 

0.5311 

6.66 

5.88 
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Table  5. 


Initial  state  of  a Freon  22  (CHCIF2)  jet  from  Figure-l-type  vessel  with  characteris- 
tics of  Table  1,  constant  initial  state  properties  of  Table  2,  and  with  rj  = 20rLj  = 
20CD^2rf^  = 0.1476m. 


Pal  = 1320kg/m^ 


P\AG  = 3.58kg/m3 
V = 232K 


t 

[S] 

^DV 

[Pa] 

V 

[m/s] 

Xj 

[kg/s] 

dMj^/dt 

[kg/s] 

0.00000 

0.4135E+07 

1.0000 

32.21 

03172 

536 

11.53 

0.00179 

03663E+07 

0.9030 

30.13 

03158 

5.01 

10.86 

0.00358 

03299E+07 

0.8117 

28.45 

03148 

4.73 

1030 

0.00537 

03008E+07 

0.7248 

27.05 

03139 

4.50 

9.83 

0.00716 

0.2769E+07 

0.6417 

25.86 

03132 

430 

9.43 

0.00895 

0.2570E+07 

0.5620 

24.83 

03126 

4.13 

9.08 

0.01073 

0.2401E+07 

0.4851 

23.92 

03121 

3.98 

8.77 

0.01252 

0.2255E+07 

0.4108 

23.12 

03117 

3.85 

8.49 

0.01431 

0.2128E+07 

03389 

22.40 

03113 

3.73 

8.24 

0.01610 

0.2016E+07 

0.2690 

21.75 

03109 

3.62 

8.02 

0.01789 

0.1917E+07 

0.2010 

21.16 

03107 

3.52 

7.81 

0.01968 

0.1829E+07 

0.1347 

20.62 

03104 

3.43 

7.62 

0.02147 

0.1749E+07 

0.0701 

20.12 

03102 

335 

7.44 

0.02326 

0.1677E+07 

0.0069 

19.67 

03099 

3.27 

7.28 
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Table  6. 


Initial  state  of  a Halon  1301  (CBrFj)  jet  from  Figure-l-type  vessel  w'ith  characteris- 
tics of  Table  1,  constant  initial  state  properties  of  Table  2,  and  with  rj  = 20rLj  = 
20CD^2rN  = 0.1476m. 


Pal  ~ 1800kg/m^ 
Pug  = 6.17kg/m3 


V = 215K 


t 

[S] 

^DV 

[Pa] 

c 

V 

[m/s] 

[kg/s] 

dMi,AL/dt 

[kg/s] 

0.00000 

0.4135E+07 

1.0000 

38.48 

0.5583 

10.73 

8.49 

0.00204 

03687E+07 

0.9030 

36.17 

0.5566 

10.08 

8.03 

0.00407 

0J338E+07 

0.8111 

34.28 

0.5552 

9.56 

7.66 

0.00611 

0J057E+07 

0.7236 

32.69 

0.5541 

9.11 

733 

0.00814 

0.2826E+07 

0.6397 

3133 

0.5532 

8.73 

7.05 

0.01018 

0.2631E+07 

0.5590 

30.15 

0.5524 

8.41 

6.81 

0.01222 

0.2465E+07 

0.4812 

29.11 

0.5518 

8.12 

6.59 

0.01425 

0.2322E+07 

0.4058 

28.19 

0.5512 

7.86 

6.40 

0.01629 

0.2197E+07 

03327 

2736 

0.5507 

7.63 

6.22 

0.01832 

0.2086E+07 

0.2616 

26.60 

0.5503 

7.42 

6.06 

0.02036 

0.1988E+07 

0.1924 

25.92 

0.5499 

7.23 

5.91 

0.02240 

0.1899E+07 

0.1248 

25.29 

0.5496 

7.05 

5.78 

0.02443 

0.1820E+07 

0.0588 

24.71 

0.5492 

6.89 

5.65 
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Figure  1.  The  experimental  arrangement. 
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p 


Figure  2.  Path  in  the  enthalpy  - pressure  plane  of  the  thermodynamic  state  of  the  initially 
liquid  agent  as  it  flows  from  inside  to  outside  the  pressure  vessel  and  achieves  a 
stable  two-phase  state  at  the  inital  section. 
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Figure  3.  Sketch  of  the  discharging  agent,  the  transition  region,  and  the  initial  section,  z 

Zj. 
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AT 


Figure  4. 


MASS  FLUX  OF  GASEOUS  AGENT 


MASS  FLUX  OF  LIQUID  AGENT 


Sketch  of  the  velocity  and  mass  fluxes  of  agent  gas  and  agent  liquid  droplets  at  the 
initial  section,  z = Zj. 
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Figure  5.  Plot  of  Table  3 results  for  Ppy  as  a function  of  time  during  vessel  discharge  of 
Freon  22. 
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Figure  6.  Plot  of  Table  3 results  for  as  a function  of  time  during  vessel  discharge  of 
Freon  22. 
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Figure  7.  Plots  of  Table  3 results  for  ^ and  Xj  as  functions  of  time  during  vessel  dischai^e  of 
Freon  22. 
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dM/dt  [kg/s] 


t[s] 


Figure  8.  Plots  of  Table  3 results  for  dMj^^dt  and  dMj^o/dt  as  functions  of  time  during 
vessel  discharge  of  Freon  22. 
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